Transport properties of 3D extended s-wave states in Fe-based superconductors 
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The Fermi surfaces of Fe-pnictide superconductors are fairly two-dimensional (2D), and it has 
thus come as a surprise that recent penetration depth and thermal conductivity measurements 
on systems of the 122 type have reported c-axis transport coefRcients at low temperatures in the 
superconducting state comparable to or even larger than that in the a6-plane. These results should 
provide important information on both the Fermi surface and the superconducting state. Here we 
consider the theory of the superfluid density and thermal conductivity in models of extended-s wave 
superconducting states expected to be appropriate for Fe-pnictide systems. We include intraband 
disorder and consider a range of different Fermi surfaces where gap nodes might exist. We show 
that recent experiments on Ba(Fei_iCoa;)2As2 can be semiquantitatively understood by such an 
approach, and discuss their implications. 

PACS numbers: 



I. INTRODUCTION 

The new Fe-based superconductors have captured the 
imagination of the theoretical and experimental super- 
conductivity community, in part because the critical tem- 
perature Tc is high in certain materials, but also be- 
cause the phenomenology of the superconducting state 
appears to be quite different from any class of novel 
superconductors heretofore discoveredii^. Since calcula- 
tions of the coupling strength Xei-ph due to the electron- 
phonon interaction show it to be small^"— , it is believed 
that electronic pairing is predominant; thus attention 
has focussed on sign-changing states with Aig symmetry 
predicted by spin fluctuation theory^"—. Depending on 
the details of the microscopic pair interaction^iHiii such 
states may display "accidental" nodal structures, nodes 
not dictated by symmetry considerations. It is important 
to understand such structures, as they provide clues to 
the origin of pairing. 

In the LaFePO, KFe2As2, and BaFe2(Asi_2;Pa;)2 sys- 
tems there is considerable evidence for low-energy ex- 
citations which could be produced by order param- 
eters with nodes, but the structure of the nodal 
manifold is still controversial. In the heavily- 
studied Ba(Fei_a;Coa;)2As2, Ba(Fei_a;Nia;)2As2 and 
BaFe2(Asi_2;Ka;)2 systems, which we will refer to col- 
lectively as Ba-122, measurements have been interpreted 
variously in terms of fully gapped behavior, deep min- 
ima in the superconducting gap, or gap nodes, and there 
is some evidence that these conclusions may be quite 
doping-dependent. The experimental probe which has 
thus far provided the lowest temperature information on 
the bulk order parameter is thermal conductivity. In the 
122 systems, the a6-plane thermal conductivity data ex- 
hibited zero or extremely small linear-T term in zero 
magnetic field, reflecting the apparent absence of any 
nodes in the superconducting gap. The field dependence, 
however, was significantly stronger than that expected 



for a large-gap superconductor—^—. This was analyzed 
by Mishra et alr^ in terms of a gap with Aig symmetry 
with no nodes but deep minima on the electron sheets. 
Bang proposed that this effect could also be explained 
phenomenologically by an isotropic Aig state with very 
small gap on one Fermi surface sheeti^. 

Recently, Martin et ali^ reported a strong linear-T 
behavior of the magnetic penetration depth for currents 
along the c-axis in overdoped Ba(Fei_i:Nia;)2As2crystals, 
compared to a predominantly behavior in the ah- 
plane, with possible much smaller ah linear-T contribu- 
tions for some dopings. These authors pointed to a need 
for theoretical analysis of superconductivity in these ma- 
terials beyond 2D models. Subsequently Reid et al^ 
measured a significant linear-T term in the low-T c-axis 
thermal conductivity of Ba(Fei_a;Coa;)2As2Crystals, com- 
pared to a smaller or zero linear term in the a6-plane. In 
addition, the anisotropy ratio {nd i^c,n) / {i^a/ K,a,N) of the 
electronic thermal conductivities normalized to their val- 
ues at Tc determined via resistivity measurements and 
the Wiedemann-Franz law, was found to increase rapidly 
as the crystal was doped away from optimal Tc. These 
authors argued that such an anisotropy ratio could arise 
only from gap nodes located precisely at flared regions of 
the quasi-cylindrical Fermi surfaces where w_F,c/w_F,a ^ 1, 
such that the corresponding arising at low T from the 
nodal structures would be negligible. 

For some samples, Reid et at reported a nonzero lim- 
iting value of Kab{T)/T of order 1 ^W/K^cm. It is im- 
portant to note that these values are much smaller than 
the value predicted in the simple BCS theory of quasi- 
particlc transport in a nodal superconductor, of order 
k/T « NoVp/{kFVA), where Nq is the Fermi level den- 
sity of states, vp is the Fermi velocity and ua is the gap 
velocity va = 9Ak/9k|nodo- This value is also much 
smaller than that observed in the cuprate case^^. The 
low-T linear term in the thermal conductivity is expected 
to be universal (disorder-independent for weak disorder) 
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in the p-wave or d-wave caso22^. Mishra et alJ^ showed 
recently that this expression continues to hold in the 
extended-s type states thought to be characteristic of 
Fe-pnictide systems, but that universality breaks down 
effectively due to the strong dependence of va on disor- 
der. 

Theory has also made some recent progress in deal- 
ing with deviations from pure 2D behavior. Early spin 
fluctuation calculations for the pairing state of the Fc- 
pnictide materials^i^ i^^'^^ found that, depending on de- 
tails of electronic structure and interaction parameters, 
nodes could occur on the electron sheets, but the gap on 
the hole sheets was always fairly large and isotropic. On 
the other hand, Graser et al^ recently presented a calcu- 
lation of spin- fluctuation pairing in the 122 systems based 
on an RPA treatment of a 5-orbital model derived from 
density functional theory (DFT) electronic structure. In 
some cases, particularly for substantial hole doping, the 
ground state was found to be of Aig symmetry, as in the 
2D case, but while the electron sheet gap was found to 
be highly anisotropic but nodeless, the hole sheet was 
found to have nodes near the top of the Brillouin zone 
near the point where the Fermi surface sheets experience 
some outward flaring. Similar results were subsequently 
obtained by Kuroki et al. for BaFe2(Asi_j,Pa;)2^. 

In this work, we consider various possible Fermi sur- 
face geometries and gap structures which may give rise to 
the unusual low transport anisotropy seen in Refs. [2TII28I . 
bearing in mind that any phenomenology which purports 
to explain the results of these works must also be con- 
sistent with the earlier thermal conductivity results of 
Refs. [iMS Although we consider 3D order parame- 
ter structures of the type found in Graser et al^, we 
do not attempt to tie our calculations or parameters to 
any particular microscopic calculation of either the Fermi 
surface or gap function, but rather to place restrictions 
on what types of structures are possible on the basis of 
these measurements. We argue that gap structures with 
small nodal segments near the = ±7r sections of flared 
Fermi surfaces are the most likely way to explain the un- 
usual anisotropy in k/T, the small size of these linear 
terms, and the magnetic field dependence observed. To 
some extent our calculations also depend on the role of 
disorder in these materials, which is not completely un- 
derstood. In Section II, we present the model we study, in 
Section III give our results, in Section IV discuss them in 
comparison with experiment, and in Section V, conclude 
and critically consider the limitations of our approach. 



II. MODEL 

A. Fermi surfaces and gap structures 

The purpose of this study is to examine transport in 
the superconducting state at very low energies; to this 
end, only the structure of the Bogliubov quasiparticle 
spectrum near gap nodes or deep minima is relevant. 



In theoretical studies of these systems performed thus 
far, nodal structures in Aig states have been obtained 
on either electron or hole pockets, but not both. We 
therefore consider two Fermi surface sheets, one of which 
{Si) possesses gap nodes. The second (52) will be as- 
sumed to have deep gap minima, or to be irrelevant al- 
together for extremely low-energy transport. Results are 
symmetrized to ensure invariance under point group rota- 
tions. For the 122 systems, we will have in mind primarily 
that sheet Si is one of the hole sheets^^ and that S2 is 
the (properly symmetrized) electron sheet. The reason 
is that, to the extent that nodes are occasionally found 
in spin fluctuation calculations on the electron sheets, 
they tend to run vertically from bottom to top of the 
Brillouin zone, and would therefore lead to an extremely 
large thermal conductivity anisotropy ^ah/'^c, in appar- 
ent contradiction to experiment. However, none of the 
results actually depend on the assignment of Si and S'2. 

Fig. [T] shows various different kinds of Fermi surfaces 
Si and nodal structures which we consider in this work. 
Among the figures labeled "^i", example 1 is a Fermi 
surface fit to the hole-doped ai sheet of Ref. [2^, calcu- 
lated using a 5-orbital fit to a first principles calculation 
for BaFe2As2 using the Quantum ESPRESSO package, 
assuming fixed experimentally determined atomic posi- 
tions in the unit cell. The nodal line structures consid- 
ered here are typical of the 3D ground states found on the 
ai and 02 sheets in that work. Example 2 is the same 
Fermi surface with horizontal nodes found in the same 
case in Ref. [2^ on the ai sheet. Example 3 is a model 
Fermi surface with identical topology but increased flar- 
ing near the top and bottom zone faces, with assumed 
V-shaped nodes near the flared portion. This Fermi sur- 
face is obtained within the same calculational scheme^^ 
for the Q!2 Fermi surface, but using a structural opti- 
mization of the internal coordinate of the As position. 
In both calculations we have used a — 3.9625 A and 
c = 13.0168 A. Similar results are found elsewhere^^r^. 
ARPES has provided some evidence for flaring of hole 
Fermi surfaces with large vp^c near the zone topt^"— , in- 
creasing with Co concentration, but at this writing there 
is no complete consensus on the 3D Fermi surface of the 
doped materials. Example 4 is the same Fermi surface 
but with horizontal nodes close to the kz value where 
fF,c takes on a maximum. Example 5 is a model surface 
similar to that found in density functional theory calcula- 
tions upon electron doping of the ai sheet together with 
a model gap structure giving nodes on the top surface. 
Finally, the last panel in Fig. [T] shows the second sheet 
S'2, fit to the Quantum ESPRESSO calculation for the 
electron sheets of BaFe2As2 for unrclaxed As coordinate, 
which is assumed to exist together with each of the cases 
1-5 for Si. We emphasize that our point of view here is 
that of phenomenologists; we allow ourselves to ask what 
would be the consequences if the Fermi surfaces actually 
took the forms shown for any of the doped materials. 

A crucial ingredient in the calculations presented below 
are the Fermi velocities near the nodal surfaces. In Table 



TABLE I: The DFT root mean square (rms) values of the in- 
plane and the out-of-plane component of the Fermi velocity 
on the different hole pockets 5*1 and on the electron pocket 82- 
In addition the rms values of the Fermi velocity components 
on the hole pockets, averaged only over the nodal regions, are 
tabulated. All values have been renormalized downward by 
a factor of 4 to account many-body effects and, are given in 
units of 10^^. In the last lines are listed the parameters 
in Eqs. [T][5] which determine the order parameter on the 5*1 
Fermi surfaces used to calculate results in this work, and the 
gap magnitude parameter Ao(meV). 



1 2 



4 5 



|T1 we list some values of Fermi velocity components and 
certain averages which will be important in the discussion 
below for the various cases for sheet Si and our fixed 
choice of 82- The dependence of the Fermi velocity 
components is exhibited in Figure [1] 

In calculations of observables presented below, we have 
reduced all Fermi velocities shown in Table 1 by a factor of 
4 to account for the effective renormalization of the bands 
with respect to DFT seen in ARPES which appears to 
be between 1.5 and 6 for all Fe-based superconducting 
materials^i^. There may be an additional low-energy 
renormalization at the meV scale, as discussed by Ben- 
fatto et alM'. All of these many body effects are lumped 
into a single doping- and momentum-independent renor- 
malization of the Fermi velocity here, to get the crud- 
est description of the anisotropy of transport properties 
driven by the underlying anisotropy of the DFT Fermi 
surface. 

The model we adopt for the order parameter is a 3D 
extension of the phenomenological form treated, e.g. in 
Ref. [1^ To obtain strictly horizontal nodes (cases 2,4,5), 
we consider gaps of the form 



A(fc_0== Ao[l + r;, cos(fc^d)]. 



(1) 



whereas V-shaped nodal structures (cases 1,3) similar to 
that found in Ref. [2^ (see Fig. [T]) are produced by 

A(fc„ ,/)) = Ao[l + r, cos(4(/.)(l - cos(fc,d))], (2) 

with > 1 for the first kind of order parameter and > 
0.5 for the second kind of order parameter. In Fig. [T]for 
simplicity we have shown only the nodal surfaces of the 
order parameters chosen on the various Fermi surfaces, 
together with the variation of the Fermi velocity on those 
surfaces as a function of k, . 
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FIG. 1: Top: various different Fermi surfaces Si cases 1-5 and 
associated gap nodal structures (dark blue lines) considered 
in this work. Top right: sheet S2 considered for all cases, with 
lines of deep gap minima (dark blue dashed lines) indicated. 

Bottom: Fermi velocity components vp^^y = y^ii|, ^ + Vp ^ 

and vp^z plotted vs. kz for the five Si cases. Fermi velocities 
used in calculations are a constant factor of 2 smaller than 
those shown here. See text for discussion. 



B. Disorder in 2-sheet pairing model 

For disorder, we will assume an orbital-independent 
matrix element which scatters quasiparticles either 
within a given band with amplitude Uu, i = 1,2, or 
between bands with amplitude Ui2. We sum all single- 
site scattering processes of arbitrary strength to obtain a 
disorder-averaged Nambu self energy E = nimpT_ , where 
nimp is the concentration of impurities, and T is the im- 
purity T-matrix as parameterized e.g. in Ref. [isl . For 
simplicity, we assume Uu = U22 = Ud, with equal den- 
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sities of states Ni = Nq throughout the paper. In our 
prehminary considerations we restrict ourselves to purely 
intraband scattering, U12 ~ 0. The disorder is charac- 
terized by two intraband scattering parameters on each 



sheet: F,- 



n-imp/i'^Ni) and 



l/{nN,Uu); For our 



simple initial case with 2 symmetric bands we set F^ = F 
and Ci = c, i = 1,2. We only consider intraband scatter- 
ing in the unitary limit^^, and further note that nonzero 
interband scattering does not affect the physics quali- 
tatively, unless the interband scattering is as strong as 
the intraband scattering^S. In presence of such strong 
isotropic scattering, one would expect a large suppres- 
sion of Tc which is not observed in experiments. 

The disorder-averaged matrix Green's function in the 
presence of scattering in the superconducting state is 
given by a diagonal matrix in band space, 



G(k,a;) 



tl'Tp -I- ekT3 -I- AkTl 



(3) 



where k = S Si is restricted to Fermi surface sheet Si 
with i ~ 1,2, and the renormalized quantities a) = w — Eq, 
Ek = Ck + S3, Ak = Ak -I- Si also depend on the band 
indices through k. The are the components of the self- 
energy proportional to the Pauli matrices Tq, in particle- 
hole (Nambu) space. 



III. RESULTS 
A. DOS 

Below wc focus on several bulk observables. It is useful 
to start with the analysis of the total density of quasi- 
particlc states (DOS) 



7VM = -^TrIm^G(k„a;), 



(4) 



where the momentum summation indicates the explicit 
integration over distinct Fermi surface sheets with mo- 
menta ki. The total residual density of states (DOS) 
at the Fermi level N{0) contains vital information about 
the low lying quasiparticle excitations and will determine 
all leading low temperature power laws. This quantity 
is sensitive to disorder or magnetic field. To give the 
reader a sense of the average density of quasiparticle ex- 
citations without reference to Fermi surface anisotropy, 
we calculate the residual DOS as a function of normal 
state scattering rate due to intraband disorder. Figure [5] 
shows A^(0) for the five different cases in consideration; 
generally it first increases with disorder, attains a max- 
imum value and rapidly drops to zero again. This "re- 
entrant" gapped behavior is a consequence of the node 
lifting phenomena characteristic of accidental nodes in an 
order parameter with Aig symmetr y ^^1'^° . Considering 
the evidence provided by Refs. [20II21I for low-energy ex- 
citations, we assume in our discussion below that we are 
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FIG. 2: Residual density of states N{iu = 0) normalized to 
the total normal state density of state at the Fermi energy 
No vs. total unitary intraband scattering rate parameter F in 
meV. Figures (a)-(e) correspond to cases 1 to 5, respectively. 



in the regime with nodes and no true gap in the system. 
The values of the scattering rate F for "dirty" systems 
presented here arc chosen in each case to correspond to 
the maximum residual DOS induced by disorder as shown 
in Fig. [21 In all cases we have verified that the corre- 
sponding Tc suppression relative to the pure system is 
0(1%) or less for aU cases. 

We note that for the various cases considered, the 
node- lifting process can take place at quite different rates. 
In particular, the nodes on the top of the closed Fermi 
sheet in case 5 are quite sensitive to a very small amount 
of disorder. The sensitivity of this case to disorder is the 
strongest argument against an explanation of the Reid et 
al. data based on a Fermi surface and state of type 5. 



B. Penetration depth 

The London penetration depth for currents flowing 
in the ith direction is a fundamental measure of the su- 
perconductivity and is related to the superfluid density 
by ps/m* = (/xoe^A^(r))~^. The superfiuid density ten- 
sor is related to the total electromagnetic response of the 
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TABLE II: Zero temperature penetration depth. 
Case Ao,ab (nm) Ao,e (nm) 



1,2 
3,4 
5 



87 
87 
87 



334 
287 
329 



system to an external field in the q 0,a; = limit. 
For a general multiband dispersion (k) and in presence 
of the impurities, one finds that the penetration depth is 
given by the expression 



1 



2Moe' 



— ( tanh — - 



dh 



(vF.i.a)^Im 



(5) 



A? 



2(vF,i.3vA,i-'5)Im 



(Af-c:;2)3/2 



where d is the distance between planes, the index i de- 
notes the band and a stands for the current direction ab 
or c and (. . . is an average over the Fermi surface. 
In the absence of disorder, this expression reduces to the 
result for the superfluid density of a clean system with 
general band found, e.g. in Ref. lU The second term 
in Eq. ([5]) contains the nodal gap slope wa = 9Ak/9k, 
which is extremely small compared to the Fermi velocity, 
and will be neglected in our calculations. Most experi- 
ments measure the change in the penetration depth AA 
from some minimum temperature Tmin, and are not sen- 
sitive to the absolute value of the penetration depth. 

Note that since we begin from a realistic Fermi surface 
and assume a Tc ~ 2bK as appropriate for the Co-doped 
122 systems, the reader might assume that we are in a 
position to calculate on an absolute length scale the zero 
temperature penetration depth Aq. Indeed the expression 
© reduces at T = to 



1 

^0 



2^06^ 

h^diT 



FS 



(6) 



which gives the correct result in the limit of a single 
parabolic band. However a closer examination of Eq. ([5]) 
reveals that the integrand determining Aq extends over all 
occupied bands, i.e. it involves the full electronic struc- 
ture and not just the Fermi velocities, whereas in Eq. (O 
we have linearized the band structure near the Fermi level 
for simplicity. Thus we are not able to use Eq. ([S]) to 
obtain accurate estimates of Aq even if the input band 
structure ei(k) was perfectly correct. By contrast, the 
integrand determining 



AA„(r) 



Ao,Q 



Ao,a 

Aq 



(7) 



is sharply peaked at the Fermi level, and so may be calcu- 
lated accurately. Even in this case, however, the overall 
scale is set by Ao,q. The usual procedure would be to 
take Aq.q from optical experiments. Indeed, in these sys- 
tems the a6-plane penetration depth is of order 200-450 
nm^"— , but there are no reported measurements of Aq.c 
of which we are aware. For purposes of this work we use 
Aq's determined from Eq. (given for completeness in 
table ini to fix the scale of AAq,. The overall length scale 
of field penetration may therefore be incorrect by up to a 
factor of two, but the T -dependence should be calculated 
accurately. 

The temperature correction AA(r) is often used as a 
probe of the nodal structure of gaps in unconventional 
superconductors^. In general, the asymptotic T — > 
power laws for systems with line nodes are AA ^ T for a 
clean and AA ^ for a dirty systen*^, but the range of 
validity of these results may be experimentally unobserv- 
able depending on the form of the impurity scattering and 
the detailed momentum dependence of the superconduct- 
ing gap. In our current calculations with intraband impu- 
rity scattering only, the concentrations have been chosen 
in order to minimize Tc suppression; the corresponding 
density of states and size of the asymptotic terms are 
therefore quite small, and the latter is not clearly visi- 
ble in Fig. [3] On the scale of the Figure shown, these 
asymptotic power laws are modified by low energy scales 
associated with both sheets Si and 5*2, as well as the sig- 
nificant energy dependence of the impurity self-energy. 
Over an experimentally relevant T range, a wide range 
of "best fit" power laws are therefore a priori possible. 

In Fig. [3l we show AA as a function of temperature 
for the various cases we consider here. The size of the 
observable quantity AA is seen to be in reasonable agree- 
ment with the measured penetration depth changes in 
recent measurements in optimal to overdoped Ba-1232^. 
The scale of A Ac is significantly larger than AAah, due to 
the much larger Ao,c- In addition. Figure |3] shows quali- 
tatively different behavior for a6-planc and c-axis pene- 
tration depths. These results, particularly those for cases 
3 and 4, are in qualitative agreement with experimental 
measurements on Ni doped 122 compounds, where lin- 
ear T behavior along the c-axis and power law behavior 
for a5-plane with exponent 1.6 — 2.5 for concentrations 
a;=0.03 to 0.07 (only those above x ^ 0.04 are outside the 
spin density wave (SDW) state)^ were reported. For Co 
doped 122 systems, wc arc not aware of any c-axis pene- 
tration depth measurements. 

The a6-planc penetration depth data in 122 systems 
has been reported for some time now to be close to a 
power law T" with an exponent a of about 2. It has 
usually been assumed that this result is obtained due 
to impurity scattering, either in a nodal systemic or in 
a fully gapped s± state with interband scattering such 
that midgap bound states are formed^. These effects 
may indeed be present in very dirty samples. However, 
the existence of a linear-T term for the c-axis penetration 
depth suggests that this simple argument cannot be the 
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FIG. 3: Change in penetration depth AA measured in nm 
vs. temperature T. Panels (a)-(e): clean penetration depth 
for c-axis (dashed lines) and afe-plane (solid lines) currents for 
cases 1-5. (f)-(j): same plots for dirty system with F chosen in 
each case to give the maximum N{u — 0) in Fig. (2] For A^t, 
power laws over the T range shown in the panel are indicated. 



complete explanation for the approximate behavior in 
the better samples; were the system sufficiently dirty to 
show behavior in the a6-plane, it would perforce man- 
ifest the same behavior for currents in the c-axis direction 
(we ignore for the moment the possibility that impurity 
scattering is itself highly anisotropic with respect to the 
ab/c directions.). The difference in power laws in the two 



directions must therefore be ascribed to some other ef- 
fect or combination of effects. In our proposal, there is no 
true behavior except perhaps at the very lowest tem- 
peratures; the upward curvature observed over a range 
of several Kelvin is due in part to the "activation" of a 
new source of quasiparticle excitations as the tempera- 
ture is increased. At low temperature all quasiparticles 
come from sheet 5*1, since we have assumed that this is 
the only sheet where true nodal excitations exist. On 
the other hand the order parameter on the 5*2 sheet is 
assumed to have deep minima A,„i„. When T exceeds 
this scale, quasiparticles from this sheet may be excited 
and a new contribution (upward curvature) to AA ap- 
pears. Since 5*2 is assumed to have smaller vp^c, it does 
not contribute significantly to A Ac, however. The second 
cause of "power law" behavior is indeed the disorder, but 
it is impossible to make universal statements about the 
large T ranges over which experiments are typically fit. 
To give a more useful measure of the form of these results 
in Fig. O we have performed best fits of the T depen- 
dence to AA ^ T" over the finite range of temperatures 
shown. 

The penetration depth anisotropy 



7A 



A^ 

Aab 



A, 



O.c 



Ao,a6 



1 + 



AA^ AA 



ab 



0,c 



A, 



0,ab 



(8) 



has been claimed in Refs. [siUs^ to be of order 3-6 at 
low temperature, and to decrease weakly with increasing 
temperature, but cannot yet be measured directly. Ex- 
perimentally, AAi is measured, and must be combined 
with measurements of Xo,i, with i = c,ab to determine 
the anisotropy ratio. Reliable measurements of Ao,c from, 
e.g. optics experiments, are not yet available for the 
122 materials to our knowledge. However, within our 
theoretical framework 7a can be determined, with the 
caveats described above, and Fig. |4] shows 7^ for the 
different cases. We therefore regard this as a prediction 
of weakly increasing "fx with temperature which may be 
verified at the point that reliable values of Ao,c are ob- 
tained from optics or elsewhere. Our results are in the 
range of anisotropy ratios deduced from experiment^, 
but increase weakly rather than decrease weakly with 
temperature at low temperatures. A glance at the ap- 
proximate inequality in Eq. ([5]) shows that a somewhat 
smaller value of Ac,o than effectively assumed in Refs. 
[Sllls^ would bring the measured AA^ 's and therefore the 
T-dependence of 7a in agreement with our results. 

We can now ask which type of Fermi surface/order 
parameter combination represents the experimental pen- 
etration depth data best. For cases 1 and 2, the c-axis 
Fermi velocity is small, which means not only that AA^ 
is smaller at low T, where it derives entirely from the Si 
nodes, but also that the contribution from the S2 deep 
minima is correspondingly more important; as seen in 
Fig. [3l for these cases the c-axis penetration depth as a 
function of T has substantial upward curvature, which is 
never observed in experiment. Power law fits to the ab- 
plane results give exponents which tend to be too close 
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presence of disorder, we evaluate the expression arising 
from the electronic heat current bubble, 



47r2reE. K(vF,i.a)Vr.) 



where 




X sech ( 



\2TJ T^' 



12 M ^ 



(m,(vFa.a)Vr,;) 



(10) 



FIG. 4: Penetration depth anisotropy 7a ^c/^ab) as a func- 
tion of temperature. Top: clean systems with F =0.0. Bot- 
tom: dirty systems with F chosen as in Fig. [S] 



to 1. On the other hand for cases 3 and 4, AAc(r) has 
no upward curvature above a scale corresponding to ap- 
proximately IK and is quasilinear over a significant T 
range, while /S.\ab{T) displays power law behavior with 
a approaching 2 in the dirty cases shown over a finite 
range, as observed in experiment. Case 5 displays a large 
penetration depth anisotropy, but an anomalous T de- 
pendence which is due to the extremely small Amm for 
this case which must be chosen to place the nodes in 
a position with high vf,c- Therefore from the point of 
view of the penetration depth temperature dependence 
and crude anisotropy, only cases 3 and 4 appear com- 
patible with experiment. We emphasize again that the 
power laws stated in Fig.[3]are phenomenological fits and 
are not exact in any asymptotic limit. Exponents may 
be expected to change as the low energy gap scales are 
changed, and in particular to increase somewhat as in- 
terband scattering is added. 

Before closing this section, we note that nonlocal elec- 
trodynamics are unlikely to play a significant role in these 
measurements. This is because the typical energy scale 
which cuts off singularities in nodal superconductors is 
^nioc = (Ca/'^a)Ao. For both directions in the 122 sam- 
ples, this scale is of order a few tenths of Kelvin or less. 
The crossover in clean samples from linear T to from 
nonlocal effects will therefore only complicate interpreta- 
tions at unobservably low temperatures, where disorder 
will probably dominate in any case. 



C. Thermal conductivity in zero magnetic field 

The electronic thermal conductivity k is calculated 
here using a straightforward extension of the standard 
approachS^i^ to multiband superconductors^. In the 



and a = a, c, is the band mass, \'F.i is the Fermi ve- 
locity on the iih band, and vf^ = Vp ^±{da^]i,i)'^ ■ Terms 
involving the gap slope va at the nodes are numerically 
negligible since va vf and are omitted in our eval- 
uations. Here, (...}</,, fc^ is an average over each Fermi 
surface sheet, denotes the sum over the bands, VF,i 
is the Fermi velocity on sheet i, and n is the number of 
FeAs planes per unit cell. 

We focus first on the linear term in k in the T — > 
limit, the only experimentally unambiguous signature of 
purely electronic heat transport by nodal excitations. In 
the hmit T 0, Eq. (P reduces to 



M 



kpVA.: 



(11) 



node 



where a is a dimensionless constant, A'nodes is the number 
of distinct nodal surfaces, assumed equivalent, and m* is 
the effective mass for motion of quasiparticles in the ab- 
plane. The result obtains exactly only if VF,i and VA,i 
are constant over the entire nodal surface, which is un- 
realistic except, crudely speaking, in a 2D situation, but 
if the velocities are replaced by the average expression of 
Eq. ([TT|) will give a reasonable result. The constant a re- 
duces to 1 /6 in the well-known case of a 2D d-wave order 
parameter, and will be of the same order of magnitude 
for a clean nodal s± state if the nodes are vertical and 
run the length of the Fermi cylinder. In any situations 
where the nodal lines do not span the dimensions of the 
Fermi sheet, i.e. consist of small circles or segments, as in 
all the cases for 5*1 considered in Fig. [1] the factor a con- 
tains the relative nodal phase space and will therefore be 
a <C 1. Note that in the d-wave and certain p-wavc cases, 
the result of Eq. ([TT|) is also universal, in the sense that 
the nodal velocities are unrenormalized by disorde r^'^i^^ , 
but in the generalized s± case this universality breaks 
dowi>i^. 

It is instructive to compare crude estimates of what 
one might expect for ImiT^o ^^i/T with actual theoreti- 
cal results and the experimentally measured values. In 
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TABLE III: Theory columns: comparison of difTerent esti- 
mates of the T ^ hmit of Eq. ^ for currents in the 
a6-plane, expressed in /iW/K^cm for a system with T'c=23K 
and d = ISA, with Fermi velocities and order parameters 
for the various cases as given in Sec. Ill Al Method Tl: 
2D d-wave universal result using {vF)rms from Table U and 
kF{vA)rnis ~ 2fcsTc aud a = 1/6, A''nodcs = 4 in Eq. PT]) : 
Method T2: Eq. (fTll using {vF,a)rms and {vA}rms averaged 
over the nodal manifold for each case. Method T3: Full 
numerical result from Eq. Ba(Fei_iCoa;)2As2: Experi- 
mental data for asymptotic low-T linear-T term in Kat from 
Ref.lUfor values of a; shown. Numbers in parentheses are val- 
ues of extrapolated linear term in /iW/K^cm for different sam- 
ples; zero is reported if these are negative. LFFO: value for 
LaFeFO from Refjp; P-122 value for BaFe2(Asi-^P:„)2with 
2;=0.33 from Ref. [53 



Theory 


Ba(Fei_ 


a:COa;)2As2 


LFPO F-122 


Case Tl T2 T3 


0.074 0.108 


0.114 0.127 


0.33 



1 118 31 1.8 

2 118 71 3.1 

3 102 31 2.5 

4 102 59 6.3 (-1,3) (-1) (-9,-13) 

5 94 47 3.1 
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Table IIIII we first shovir progressively more accurate es- 
timates of this quantity. First, wc remark that a small 
lincar-T term in Kat must be present if nodes are present 
and if Vi? has any nonzero component in the a6-plane. 
To be consistent with earlier data and with the recent 
work of Reid et ali^ in the Co-doped Ba-122 system, this 
term must be quite small, consistent with zero within 
the resolution of the experiment, at least near optimal 
doping. It is clear from the table that a 2D d-wave 
state or, equivalcntly, an anisotropic nodal s-wave state 
with line nodes running the length of the Fermi cylin- 
der ("method Tl") would dramatically overestimate the 
experimental results for k/T in the ab plane. The best es- 
timate based on Eq. pT|) ("method T2"), accounting for 
realistic values of nodal velocities wa, still gives a value 
significantly larger than experiment near optimal dop- 
ing in the Ba(Fei_2;Coa;)2As2 system studied by Reid et 
al. as shown in Table [1111 This is at least in part because 
method T2 also does not account for the small nodal 
phase space. The T3 column of Table [Till shows the clean 
limit of the exact numerical evaluation of Eq. ^ for the 
5 cases. These results are now of order the scatter in the 
data (consistent with zero) of the a6-plane linear-T terms 
in K identified in samples near optimal doping by Reid et 
al., which we take to be the resolution of the experiment, 
about ^ l^W/K^cm. 

Depending on the Fermi surface, the corresponding c- 
axis transport can be smaller or larger than the a5-plane 
results discussed. To illustrate the qualitative anisotropy 
found here, we plot in Fig. [5l the a6-plane and c-axis 
thermal conductivities k/T for increasing intraband im- 
purity scattering. There are several interesting points in 
this Figure. First, we note that for cases 1 and 2, where 
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FIG. 5: Residual thermal conductivity limT^o'v(r)/r in 
units of /i W/K^cm as a function of impurity scattering rate 
F in meV. Panels (a,c,e,g,i) correspond to absolute thermal 
conductivity for thermal current in the a&-plane (solid) and 
along the c-axis (dashed), respectively, for cases 1-5 depicted 
in Fig. [T] Panels (b,d,f,h,j) correspond to normalized thermal 
conductivity k/kat for the same cases. 



the Fermi surfaces are nearly cylindrical, the a6-plane 
conductivity is much larger than the c-axis conductivity. 
Clearly again these cases are inconsistent with experi- 
ment. For the more flared Fermi surface, cases 3 and 4, 
it is possible to have Kab/T and Kc/T comparable, or even 
have Kc > Kab- Note this is more likely for the horizontal 
line nodes of case 4, which have been chosen to lie at the 
maximum value oi vpc for this Fermi surface, whereas 
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FIG. 6: Normalized thermal conductivity anisotropy ratio 
7k = \vaiT^o{i^c / i^c,n) / {nab / K-ab,N) as a function of the disor- 
der scattering rate normalized to the critical value of disorder 
Fnl that lifts the nodes. 

the V-shaped nodes of case 3 sample the high vp^c values 
less often. A large c/a anisotropy is also observed in case 
5 for reasons similar to those of case 4. 

As intraband scattering is increased, a constant term 
is added to the off-diagonal self-energy in Nambu spaco^^, 
which eventually lifts the nodes, as seen in the DOS 
(Fig. [2]). As this occurs, two effects happen in parallel: a) 
wa decreases; and b) the nodal surfaces (which are shown 
in Fig. [1] for the clean system only) move towards the top 
and bottom of the Brillouin zone, respectively. The lat- 
ter effect can increase or decrease the thermal conductiv- 
ity depending on whether or not the node-lifting process 
moves the nodes towards regions of higher or lower vf,c- 
Thus, depending on details, the c/a anisotropy can be 
enhanced or suppressed by disorder. 

In Ref. [21I, the authors specifically discuss the normal- 
ized measure of thermal conductivity anisotropy 7^ = 
limT->.o('«c/ Kc,N)/{Kab/ Kab,N), and demonstrate that this 
ratio has a minimum at optimal doping. This result 
should be treated with caution, however, both because 
a) the linear term in Kat is consistent with zero for much 
of the doping range; and b) the normal state kjv's are de- 
termined by measurement of the resistivity extrapolated 
to T = and use of the Wiedemann- Franz law. However, 
if the elastic transport scattering rate 1 /rtr is reasonably 
isotropic, as seems likely, 7k will be determined simply 
by the ratio 

((4,c)node/(4e))/((4afc)nodc/(4,ah», (12) 

where {vp Q,)nodc is the average of the Fermi velocity com- 
ponents a over the nodal manifold. Here we calculate 
7k within a simple model where normal state heat con- 
duction is controlled by (isotropic) intraband impurity 
scattering only. In Fig. [SJ we show that only for cases 
3-5 it is possible to achieve values of 7^ substantially 
greater than 1. For case 5, the node lifting scattering 



rate Fnl is minuscule in absolute units, such that we 
rule out this case as well, provided our assumption that 
intraband scattering dominates is correct. 

An alternate possibility to create a small residual k/T 
would be the creation of an impurity band in a gapped 
state due to interband scatteringi^, but in this case there 
is no natural reason to have a transport anisotropy other 
than that given by Kc/f^ab ~ {"^f c) / ("^f ab) ^ which for all 
cases discussed here would give a very strong dominance 
by the a&-plane thermal currents on an absolute scale, 
and would lead to 7^ ~ 1 in all samples, in contrast to 
experiment. 

D. Magnetic field dependence of k, 

The final notable feature seen in experimcn t^^i^^ is an 
unexpectedly large enhancement of the in-plane residual 
k/T term in the presence of a magnetic field along the 
c-axis. With moderate magnetic field of a few Tesla, the 
residual in-plane term K.ab{H,T 0)/T becomes com- 
parable to the c-axis term Kc{H,T 0)/r, and both 
exhibit a field dependence with downward curvature rem- 
iniscent of the presence of field-induced residual quasipar- 
ticle excitations contributing to the thermal current^. 
This effect was discussed already for the a&-plane ther- 
mal conductivity in Ref. [isl where it was concluded that 
the most likely origin of the strong field dependence in 
apparently fully gapped samples was contributions from 
deep minima in the gap on the electron sheets, as found 
in early multiorbital spin fluctuation theories^ i^^'^^i^^'^^ . 
Note that in the 3D spin fluctuation analysis of Graser et 
ali^, these deep minima on the electron sheets continue 
to exist, and run the entire length of the Fermi cylinders, 
thus providing much more phase space than at zero fleld, 
as soon as the magnetic fleld energy spans the very small 
minimum gap Amin- Thus the ab thermal conductivity, 
negligible in zero field, can jump fairly rapidly^. 

To include the presence of magnetic field explicitly, 
in this work we use the standard semiclassical (Doppler 
shift) metho d^^'^"" — , where the quasiparticle energy is 
modified by its Doppler shift in the local superfluid flow 
fleld. While this approach is never exact^, it provides 
a good qualitative description of the magnetic fleld de- 
pendence of the quasiparticle properties of nodal super- 
conductors over a wide fleld rango^. The local Doppler 
shift of a quasiparticle with momentum k (band indices 
are suppressed) is 

,5c^(k,r) =vj.(k)-p,(r) (13) 

where (r) is the local supercurrent momentum, approx- 
imated for a single vortex as 

Ps(r)-f^ (14) 

Here r is the distance from the vortex core and 9 is the 
vortex winding angle. In the presence of the magnetic 
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FIG. 7: Residual thermal conductivity k/T in units of 
/iW/K^cm vs. magnetic field H in Tesla. Panels (a)-(e): 
cases 1-5 for Fermi surface 5*1 with intraband scattering rate 
r = 1.3 meV (Case 1), F = 0.6 meV (Case 2) and F = 1.9 
meV (cases 3,4), and F = 0.06 meV for case 5. 



field, we replace w by w — (5a;(k, r) in Eq. In this ap- 
proach, any physical quantity is calculated in the usual 
way with Dopplcr shifted energy, and the averaging is 
performed over a unit cell of the vortex lattice, approxi- 
mated as circular, as 



k{H) 



1 



Rh r^TT 

drr / deK{r), (15) 



where K(r, 6) is the local (Doppler shifted) result from 
Eq. Rh — v^o/tt/? is the magnetic length scale, 
with $0 the superconducting flux quantum. For lower 



fields H <C i?c2, this is a qualitatively good approxi- 
mation. Somewhat more sophisticated approximations 
can be obtained within the semiclassical framework^, 
but the above will suffice for our qualitative purposes. 
Here Hc2 = $o/(27r^Q) is the upper critical field and 
^0 = vf/{'2.'kTc) is the coherence length. 

Using this method, we now calculate the low-T resid- 
ual thermal conductivity in both directions as a function 
of the magnetic field. Fig. [7] shows the field dependence 
of the residual linear term, which is similar to a supercon- 
ductor with nodes, for both directions of the heat current, 
with the magnetic field taken always to point along the 
c-axis as in experiment. Note that in principle the sit- 
uations with heat current parallel and perpendicular to 
the field are physically different regardless of the orien- 
tation of the crystal axes: in one case the heat current 
fiows down the axis of a vortex and may be assumed to 
sample a single ^(r), while in the other case it samples a 
modulated k{y) across the vortex lattice. This averaging 
problem for the perpendicular heat flow has never been 
addressed in a completely convincing fashion. We there- 
fore present two extremes of the possible fleld dependence 
for currents in the the a&-plane: first, a "parallel average" 
over the vortex flow field (/tah) = (^(r)), where (...) rep- 
resents a direct average over the vortex unit cell as in 
Eq. (fTS]): second, a "series average", Rab = ('^(r)"^)"^, 
as discussed in Ref. [H. We note that the results for the 
Brandt- Pesch-Tewordt average Green's function theory^ 
appear to lie closer to the parallel average. 

Fig. [7] shows several interesting effects. First of all, 
while the Fermi surfaces with less flaring (cases 1 and 
2) displayed much larger values of the absolute a6-plane 
thermal conductivity in zero fleld (Fig. [SJ panels (a) and 

(b) ) than the more flared cases 3 and 4 (Fig. [SJ panels 

(c) and (d)), the afe-plane field dependence Ka{H) at any 
reasonable field is seen to dominate the absolute c-axis re- 
sponse Kc{H) for all cases 1-4 (Fig. [71 panels (a-d)). This 
is because we have assumed, as in Ref. [39l that the ^2 
sheet has a deep gap minimum, spanned by the Doppler 
field energy Eh at a field H of only 1.2 T. Note this is 
an average scale over the entire sheet; the dispersion of 
the Fermi velocity with leads to a smearing of any 
feature associated with this scale. It is these long deep 
gap minima on 5*2 (putative electron sheets) which dom- 
inate the field response in this case, with subdominant 
contributions from the 5*1 (putative hole sheets) above 
this field range. Because Kcn ^ i^aN, it is still possible 
to have a normalized thermal conductivity larger in the 
c-axis direction (Fig. [3 panels (f) and (g)). 

The Rcid et al. experiment Ref. I2l| shows a pronounced 
tendency for the normalized anisotropy ratio 7„ to fall 
rapidly from large values as a function of magnetic field, 
and then to saturate at values of order 1 when fields of a 
few Tesla are reached. Within our calculational scheme, 
such a conclusion is quite natural since the zero-field 
anisotropy (Fig. [B]) 7/^ is dominated by the weak nodes 
on the Si sheet(s), which can give large 7k's, particularly 
in cases 3 and 4, as discussed above. On the other hand. 
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FIG. 8: Inverse normalized thermal conductivity anisotropy 
ratio = limT->o('^a/'^a,jv)/(Kc/Kc,jv) a function of mag- 
netic field for Case 3. Solid line: parallel averaging (nab)', 
dashed line: series averaging ((k^^^))"^ as described in text. 
Symbols: data for highly overdoped x = 0.127 sam- 

ple of Ref. l2ll . which shows nonzero Kab/T as T — > 0. 



as the field is switched on, the contributions of the deep 
gap minima on the ^2 sheets come into play, and indeed 
dominate for sufficiently large fields because of the much 
larger phase space of the gap minima on iS'2 compared 
to the phase space taken by the gap nodes on ^i. If a 
single 2D sheet dominates the field dependence, the nor- 
malized anisotropy ratio will reduce to approximately 1 
since the Fermi velocity anisotropy has been removed by 
the normalization, provided the averaging over the vor- 
tex lattice is performed identically for in-plane and c-axis 
currents. This behavior is seen in FigurelH where the in- 
verse of the anisotropy ratio, 7^^ is plotted to facilitate 
comparison with experiment. If scries- type averaging is 
more appropriate in the a6-plane, the high-field 7"^ can 
be considerably larger, as also seen in Figure [5] 



2.2 2.4 26„ 



3 3.2 3.4 3.6 



FIG. 9: Normalized thermal conductivity k/kjv for heat cur- 
rent in a and c directions vs. degree of flaring of Fermi surface 
Si, parametrized by ^;™'"'/'^;^"^ 



tivity as outlined above. The superconducting order pa- 
rameter is held fixed in this process. Within our current 
interpretation of the data on Co-doped Ba-122, increas- 
ing the amount of fiaring is equivalent to doping the sys- 
tem, although the correspondence between the two scales 
is nonlinear and not completely determined. Neverthe- 
less, it is clear that an essential feature of the experiment, 
the strong increase of the normalized anisotropy ratio 7k 
can be understood by making this simple assumption. 
Note that the normalization automatically divides out 
the average effect of the increasing VF,c/vF,ab] what hap- 
pens, therefore, is that the flaring increases particularly 
near the nodes and therefore preferentially enhances low- 
energy c-axis quasiparticle transport. Fig. [9] illustrates 
the effect of the increased flaring, parameterized by the 
ratio of the Fermi velocity components, for the order pa- 
rameter with V-shaped nodes as in cases 1,3. The effect 
is superficially quite similar to the effect of doping on the 
anisotropy ratio observed by Reid et al^. 



E. Doping dependence 

Within our current crude approach, doping is accom- 
plished by rigidly shifting the chemical potential rela- 
tive to the parent compound. However, the dependence 
of the Fermi surface thus obtained is consistent neither 
with that observed in ARPES experiments on 122 com- 
pounds nor with virtual crystal approximation within 
DFT calculations. In particular, the flaring of the hole 
(Si) Fermi surface upon electron doping docs not in- 
crease rapidly with electron doping, as seen in ARPES 
and DFT calculations^!^. In this section we therefore 
investigate, independently of any microscopic model of 
the band structure, the effect on superconducting state 
transport of the flaring of the Si Fermi surface. To this 
end, we simply interpolate the Fermi surfaces between 
cases (1,2) and (3,4), and calculate the thermal conduc- 



F. Modulations of k with respect to magnetic field 

angle 

Measurements of the speciflc heat and thermal con- 
ductivity as a function of the angle between an applied 
magnetic field and the crystal axes are in principle power- 
ful methods to determine the location of gap node o^^'^^ . 
Since these measurements are probes of the bulk prop- 
erties, and the treatment of surfaces in the Fe-pnictide 
systems continue to pose problems, rotating field mea- 
surements were proposed early on as appropriate ways 
to map nodal structures of the order parameter in these 
systems^. Measurements of the low temperature spe- 
cific heat as a function of the angle of an in-plane mag- 
netic field have been reported very recentlj*^. While this 
method is sensitive to vertical line nodes, it misses hor- 
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FIG. 10: Residual specific heat C(H)/C]v (top) and c-axis 
thermal conductivity (bottom) Kc(H, T — >■ 0)/kcn normalized 
to zero field normal state values as a function of direction of 
applied magnetic magnetic field described by polar angle 6z 
and azimuthal angle $. The Fermi surface and form of order 
parameter used for this calculation is same as case 3 with 
Ah = 0.3meV and Ae = l.SmeV. Magnetic field strength is 

lo.or. 



izontal nodes altogether and can also miss nodal struc- 
tures near the extrema of the Brillouin zone for a flared 
Fermi surface. Similar considerations apply to the oscil- 
lations of the thermal conductivity with field angle, which 
historically have been more sensitive to weak nodal struc- 
tures. Here we have calculated, for an arbitary direction 
of the magnetic field, the thermal conductivity for a fixed 
amplitude of the magnetic field for case 3. 

The results of these calculations are shown in Fig. (TU] 
The situation for polar angles 9z near zero, i.e. for an 
alignment of the field along the c-axis direction and the 
heat current applied with an in-plane angle $, clearly 
causes negligible thermal conductivity oscillation. As the 
polar angle is increased, the expected characteristic four- 



fold oscillations emerge^. Since we have evaluated the 
conductivity at T — >■ 0, minima as a function of $ for con- 
stant 9z reflect the existence of nodes at angles 90° — 9z- 
We have evaluated the T dependence as well, and find 
that inversions of this pattern can occur at higher T—'^—, 
but have not focussed on this aspect here. Within our 
model, for most angles the observed fourfold oscillations 
with minima at 0° and 90°, correspond to the deep min- 
ima of the gap function on the 5*2 sheet assumed in our 
model. Only when 9z approaches 90° are the nodal sur- 
faces on Si probed by the field. One should expect for 
cases 1 and 3 (V-shaped nodes) new double minima on 
both sides of $ =0°, 90°, etc. to appear due to the 
V-shaped nodes for polar angles 9z ^ 90°. These are 
not visible in the Figure, however, in part because the 
anisotropy of the Fermi surface masks small nodal fea- 
tures. In our model it appears difficult to observe these 
small features, but it may nevertheless be possible to de- 
tect them with high experimental resolution. Since the 
cases with full horizontal nodes are robust, one might 
think that a polar angle scan of the magnetic field would 
easily pick up minima due to these features, but this 
effect, while present, is obscured in our calculations by 
the c-axis dispersion of the Fermi velocity. Our study 
suggests that a careful mapping of the weak Si nodal 
structures with this method, while clearly desirable, may 
unfortunately be difficult both due to the small size of 
the related signals and of course the experimental imple- 
mentation of the 3D field rotation. 



IV. CONCLUSIONS 

In the presence of an order parameter with 3D struc- 
ture on an anisotropic FS, transport properties may be 
qualitatively very different along the plane and the c axis, 
and FS topology may also play an important role. In 
particular, the Fermi velocity is a characteristic weight 
factor for transport coefficients: a larger Fermi velocity 
near the high density of quasiparticles around gap nodes 
will greatly enhance the superfluid density, microwave 
or thermal conductivity. Our calculations, within the 
framework of anisotropic s± states found in recent spin 
fluctuation calculations^^ for 122 systems, support the 
suggestion of Reid et al. that the anisotropy observed in 
the thermal conductivity at low T in 122 systems is most 
likely due to the occurrence of gap nodes on portions 
of the Fermi surface with large c-axis component of the 
Fermi velocity. It is noteworthy that from the point of 
view of microscopic theory, the nodal structures which 
dominate the c-axis transport in the Co-doped Ba-122 
materials are most likely to be found on the hole rather 
than electron Fermi surfaces, in contrast to calculations 
appropriate for 1111 and other more 2D materials which 
found the opposite. It is also possible that changes in mi- 
croscopic interaction parameters can lower the gap min- 
ima on the electron sheets so as to become nodes, as may 
occur in the P-doped 122 material. Our calculations are 
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phcnomenological in nature and have only assumed that 
the nodes occur on one sheet of either type. 

Linear-T terms in the thermal conductivity can arise 
whenever a finite density of residual quasiparticle states 
exists at the Fermi level due to a small amount of dis- 
order. Some suggestions have been made that experi- 
ments indicating low-energy quasiparticle states can be 
understood by considering an isotropic s± state without 
nodes but with interband scattering, such that an impu- 
rity band is created at the Fermi level^. At least in the 
122 systems, this possibility now appears to be conclu- 
sively ruled out, for several reasons. First, such an im- 
purity band would lead to a dominant dependence of 
the penetration depth, not linear-T, as reported by Ref. 
[20I in 122 systems). Secondly, at low T disorder scatter- 
ing is the dominant scattering mechanism and also de- 
termines the characteristic transport time. Calculations 
with realistic impurities show weak variation of the re- 
laxation rate around the Fermi surface; and the density 
of impurity band states is therefore roughly isotropic; 
thus transport anisotropy will be determined in such a 
situation entirely by the ratio {vp^)/{vp^f) <?C 1, rather 
than by the Fermi velocity values on special parts of the 
Fermi surface. This result is now inconsistent with Ref. 
[21I . For this reason we have neglected a treatment of 
interband impurity scattering in this work, but realistic 
calculations will require a knowledge of both intra- and 
interband scattering amplitudes of the dopant impurities. 

In our work here we have also neglected magnetic im- 
purities, in part because the microscopic calculations of 
Kemper et. ali^ reported a weak magnetic potential 
component of the Co dopant. In addition, magnetic im- 
purities are generically strong pairbreakers in Aig states. 
Similar to the discussion of interband nonmagnetic scat- 
terers above, if they were strong enough to form impurity 
bands on the Fermi level and thus influence transport, 
they would create significant Tc suppression, large val- 
ues of the linear residual term in Kat, and a normalized 
anisotropy 7^ in zero field close to unity, none of which 
are observed in experiment. 

We have focussed here primarily on two properties, 
penetration depth and residual thermal conductivity in 
the T ^ limit. For two concrete examples of Fermi 
surfaces "inspired" by forms which occur in density func- 
tional theory calculations for certain 122 systems, we 
have shown that one can simultaneously explain the exis- 
tence of large penetration depth and thermal conductiv- 
ity anisotropics 7a , 7k > 1 , indicating a predominance of 
nodal quasiparticles with large c-axis components of the 
Fermi velocity. Because of the complicated multiband 
nature of the system, as well as because the order pa- 
rameters assumed have important low-energy minimum 
gap scales, asymptotic low T power laws will rarely be 
observed in the penetration depth AA(r) if the T de- 
pendence is fit to an experimentally relevant range. We 
have therefore tried to give fits AA(r) ^ T" to our the- 
oretical results over such a finite range, and found for 
two promising cases with "flared" or closed Fermi sur- 



faces in-plane exponents aab ^1.5 — 2 which depend on 
details of order parameter structure and on topology of 
Fermi surface. On the other hand, the c-axis penetra- 
tion depth change is very close to linear T dependence 
for these cases. These results are quite close to those 
obtained in experiment by Martin et al^. Quantitative 
comparisons appear to us to be impossible at the present 
time, however, due to uncertainties in the magnitudes of 
the zero T penetration depths. 

In the case of thermal conductivity, only the coeffi- 
cient of the linear T term is relevant since it can be un- 
ambiguously assigned to electronic transport. We flnd 
again in the two promising cases that qualitative thermal 
conductivity anisotropy in agreement with experiment is 
found. For this calculation, semiquantitative agreement 
with experiment can be obtained once the disorder scat- 
tering rate and overall Fermi velocity renormalization is 
flxed. One further remarkable aspect of the experiments 
explained by the 3D Fermi surface and order parameter 
models we have considered is that the linear-T ab-plane 
thermal conductivity is substantially smaller than one 
might expect on the basis of simple dimensional analy- 
sis, reflecting the small phase space from the short nodal 
segments near the top of the Brillouin zone. The c-axis 
measurements of Reid et ali^, in this picture, force one 
to the conclusion that Kab is not zero in these materials, 
as previously reported, but must be so small as to be 
unrcsolvable within current experiments. 

Further conclusions can be drawn from study of the 
magnetic field dependence of the residual linear terms 
within the semiclassical approach. The observed field 
dependence in the 122 is quite similar to nodal systems 
in clean limit, including the size of the term. This sug- 
gests that the small nodal structures probed by the c-axis 
transport are not giving rise to the field dependence; in- 
stead, this contribution is dominated by the remaining 
Fermi surface sheets which have no nodes but rather deep 
gap minima which run the length of the Fermi cylinder—. 
An immediate consequence is that in the cases which ap- 
pear promising relative to experiment in zero field, the 
a6-plane transport is immediately enhanced upon appli- 
cation of a small field, leading to comparable nodal-type 
field variation of the thermal conductivity in both direc- 
tions, as reported in Ref. [2ll . 

Finally, we have briefiy discussed how the nodal struc- 
tures which appear to us to be the most likely in this 
3D system could be detected. We propose a full angle- 
dependent measurement of the directional thermal con- 
ductivity. Weak 4-fold oscillations of the specific heat 
have been reported in other Fc-bascd superconductors^, 
and discussed theoretically rcccntlyiiS On the other hand, 
if the heat current is along c and the field is always in 
the plane, the unique features of the nodal structures 
near the top of the Brillouin zone face will be missed. 
Therefore measurements at polar angles of the magnetic 
field between and 90 degrees with respect to the plane 
may be useful. Our study has found, however, that the 
angular dependence of this signal is probably dominated 
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by the structures on electron, not hole sheets; therefore 
these small hole pocket signatures may be difBcult to ob- 
serve. 

Our calculations provide further evidence for the con- 
jecture of an "intrinsic sensitivity" of the gap structures 
in the Fe-based superconductors to small perturbations 
in electronic structureii. Recall that according to Graser 
et al^, the nodes which appear near the top of the 
Brillouin zone on the hole pockets in these systems are 
probably related to the suppression of d^z/dyz orbital 
weight there, which in turn occurs because of special fea- 
tures in the 122 crystal structure not shared by other 
Fe-pnictides. 

Within our model, the increase of normalized thermal 
conductivity anisotropy observed by Rcid et al—- in Co- 
doped Ba-122 may be attributed to the increased flaring 
of the Fermi surface sheet as seen in ARPES experiments 
and DFT calculations. This is a nontrivial effect due 
not to the increase of the average c-axis Fermi velocity 



but the increase of this component at the position of the 
gap nodes. There are clearly other effects present, how- 
ever, including the possible occurrence of nodes on the 
electron sheets for the most overdoped systems, as indi- 
cated by the non-zero Kab- It is an important challenge 
to construct a microscopic theory which can simultane- 
ously describe the evolution of the Fermi surface and gap 
structure with doping. 
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